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Abstract 

We provide a Hamiltonian analysis of the Mixmaster Universe dy- 
namics on the base of a standard Arnowitt-Deser-Misncr Hamiltonian 
approach, showing the covariant nature of its chaotic behaviour with 
respect to the choice of any time variable, from the point of view either 
of the dynamical systems theory, either of the statistical mechanics one. 

1 Introduction 

Since Belinski-Kalatnikov-Lifshitz (BKL) derived the oscillatory regime which 
characterizes the behaviour of the Bianchi type VIII and IX cosmological 
models Q (the so-called Mixmaster universe Q) near a physical singularity, 
a wide literature faced over the years this subject in order to provide the 
best possible understanding of the resulting chaotic dynamics. 

The research activity developed overall in two different, but related, di- 
rections: on one hand the dynamical analysis was devoted to remove the 
limits of the BKL approach due to its discrete nature, on the other one to 
get a better characterization of the Mixmaster chaos (especially in view of 
its properties of covariance). 

The first line of investigation provided satisfactory representations of the 
Mixmaster dynamics in terms of continuous variables (leading to the con- 
struction of an invariant measure for the system || , [|| ) . 
The efforts in the second direction found non-trivial difficulties due to the 
impossibility, to apply the standard chaos indicators to relativistic systems. 

The existence of difficulties related to the covariance of various ap- 
proaches prevented, up to now, to say a definitive word about the Mixmaster 
chaoticity (apart from the indication provided by 0], see also @). 

The aim of this work is to show how in the points of view of the theory 
of dynamical systems and statistical mechanics the representation of the 
Mixmaster chaoticity is independent on the choice of a time gauge. 



The description of the system evolution as a "stochastic scattering" iso- 
morphic to a billiard on the Lobachevsky space can be constructed inde- 
pendently of the choice of a time variable, simply providing very general 
Misner-Chitre-like coordinates. 

On the other hand, we show how the derivation of an invariant measure 
for the Mixmaster model (performed in || ^J] within the framework of the 
statistical mechanics) can be extended to a generic time gauge. In fact, 
asymptotically close to the cosmological singularity, the Mixmaster dynam- 
ics can be modeled by a two-dimensional point-Universe randomizing in a 
closed domain with fixed "energy" (just the ADM kinetic energy); since it is 
natural to represent such a system by a microcanonical ensemble, then the 
corresponding invariant measure is induced by the Liouville one. 

Up to the limit of the adopted approximation on the form of the poten- 
tial term, our analysis shows, without any ambiguity, that the Mixmaster 
stochasticity can not be removed by any redefinition of the time variable. 

2 The Hamiltonian Formulation 

The geometrical structure of the Bianchi type VIII and IX spacetimes, i.e. of 
the so-called Mixmaster Universe models, is summarized by the line element 

i 

ds 2 = -N{r,fdrf + e 2a a { a j (1) 

where N(r]) denotes the lapse function, a' 1 are the dual 1-forms associated 
with the anholonomic basis [] and fyj is a traceless 3x3 symmetric matrix 
diag(/?n, /?22, P33)', a, N, fiij are functions of r] only. Parameterizing the 
matrix (5ij by the usual Misner var iables § 

0u = 0+ + V3/3- , P22 = (3+ - Vs/3. , (3 33 = -2f3+ (2) 

The dual 1-forms of the considered models are given by: 
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the dynamics of the Mixmaster model is described by a canonical variational 
principle 51 = 5 J L dij = 0, with Lagrangian L 
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Here ()' = D = dete a+/3 ^ = e 3a and the potential V (a, /?+,/?_) reads 



V = - (.D 4 ^ 1 + D iH2 + D m ' A ^j - D 2Hl+2H ' 2 ± D 2H2+2Hs ± D 2Hz+2Hl , (4) 

where (+) and (— ) refers respectively to Bianchi type VIII and IX, and the 
anisotropy parameters Hi (i = 1,2,3) denote the functions || 

1 /3 + + V3/?_ 1 0+-V3P- 1 2/3+ 

H\ = — , -Ho = — i = • (5) 

3 3a 3 3a 3 3a 

In the limit Z? — > the second three terms of the above potential turn 
out to be negligible with respect to the first one. Let's introduce the new 
(Misner-Chitre-like) variables 
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with / denoting a generic functional form of r, 1 < £ < oo and < 9 < 2-7T. 
Then the Lagrangian (||) reads 



TV 



Zv(f{T),t,0). (7) 



In terms of / (r), £ and we have 

D = ea;p{-3£- e /(r) } 



(8) 



and since D — > toward the singularity, independently of its particular 
form, in this limit / must approach infinity. The Lagrangian (||) leads to 
the conjugate momenta 
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which by a Legendre transformation make the initial variational principle 
assume the Hamiltonian form 



S J (ft? +PeO' +p T 



Ne- 2f \ 

T -^p- H )<*n = o, (io) 



being 



H —^4 + « 2 (e - l) + ^ + 24 Ve". (11) 



(*)' 



By variating (fTO) with respect to iV we get the constraint H = 0, which 
solved provides 

-Pr = ^- Hadm = ^ • v^ 2 + 24Fe 2 / (12) 
ar ar 

where 

e a = (e a -i)pe 2 + ^ I . d3) 

In terms of ((l^) the variational principle (|i"o| ) reduces to 

^ / (p^ + P^' - f"H A DM) d V = 0. (14) 

Since the equation for the temporal gauge actually reads 

12D ofdf , 
it ADM dT 

our analysis remains fully independent of the choice of the time variable 
until the form of / and r' is not fixed. 

The variational principle ( |i~4| ) provides the Hamiltonian equations for £' 
and & | 

'-sL^- 1 )*. '=«L«£iT <16) 



Furthermore can be straightforward derived the important relation 
d(H A DMf') _ djHADMf') _^ djHADMf') _ d(H A DMf') 

dr] drj df df 

i.e. explicitly 



(17) 



&HADM ^ -24.f2F + ^l. (18) 



df ^Hadm V df 



2 In this study the corresponding equations for p'^ and p'g are not relevant. 
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In this reduced Hamiltonian formulation, the functional / (77) plays simply 
the role of a parametric function of time and actually the anisotropy param- 
eters Hi (i = 1, 2, 3) are functions of the variables £, 9 only 



U = 1 - ^ - fcos0 + V3sin( 
3 3£ V 



3 3£ 

H 2 = 1- (cos 0-^/3 sing) (19) 

tf 3 = i + 2^HIcos.. 

Finally, toward the singularity (D — ► i.e. / — > 00) by the expressions 
(gg HI), we see that 

|j = O (efv) . (20) 

Since in the domain Tjj ah the Hi are simultaneously greater than 0, the 
potential term U = e 2 -^V can be modeled by the potential walls 



U OB = 8^ *)) + ©00 + (21) 

+00 if x < 
i/ x > 



6oo (a?) 



therefore in Th the ADM Hamiltonian becomes (asymptotically) an integral 
of motion 

ft\aT \ HaDM = + 2A ' U = £ = E = const 

\/{£,V\£L H < an adm _ dE _ n ( 22 ) 

I df ~ df ~ U • 

The key point for the use of the Misner-Chitre-like variables relies on 
the independence of the time variable for the anisotropy parameters Hi. 



3 The Jacobi Metric and the Billiard Representa- 
tion 

Since above we have shown that asymptotically to the singularity (/ — > 00, 
i.e. a — ► —00) cLTLadm / df = i.e. TLadm = e = E = const., the variational 

3 By 0() we mean terms of the same order of the inclosed ones. 
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principle ( [L4D reduces to 

5 J (pfdC + p e d# - Edf) = 5 J (p € d£ + p*d0) = , (23) 

where we dropped the third term in the left hand side since it behaves as 
an exact differential. 

By following the standard Jacobi procedure 1C ] to reduce our variational 
principle to a geodesic one, we set x al = g ab Pb, and by the Hamiltonian 
equation we obtain the metric 

/e-^ (<"-!). 9 >> - , ~- P4) 

By these and by the fundamental constraint relation 

(e-l)n 2 + ^- = E\ (25) 



we get 

9ab x a 'x b ' = t I (e 2 - l) Pi ? + j£j ) = f'E. (26) 



E v 

By the definition x a ' = £j = u a f- , the @) rewrites 



= /'iS, (27) 



which leads to the key relation 



d^S^*. (28 ) 

Indeed the expression (|2S|) together with + = E 1 /' allows us to put 
the variational principle (23) in the geodesic form: 

5 J f'E dr, = 6 J ^g ab u a u b f'E ds = 5 J ^G ab u a u b ds = (29) 

where the metric G ab = f'Eg ab satisfies the normalization condition G ab u a u b = 
and therefore ^ 

ds „ „. (is _ , . 

4 We take the positive root since we require that the curvilinear coordinate s increases 
monotonically with increasing value of /, i.e. approaching the initial cosmological singu- 
larity. 
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Summarizing, in the region Th the considered dynamical problem reduces 
to a geodesic flow on a two dimensional Riemanniann manifold described by 
the line element 



ds 2 



E 2 



+ 



(« 2 



i )de 2 



(31) 



The above metric (31) has negative curvature, since the associated curvature 
scalar reads R = — -Jj; therefore the point-universe moves over a negatively 
curved bidimensional space on which the potential wall (Q) cuts the region 
Th- By a way completely independent of the temporal gauge we provided 
a satisfactory representation of the system as isomorphic to a billiard on a 
Lobachevsky plane flOfl , 

The invariant Lyapunov exponent defined as |11] reads 



A„ = sup lim 



In Z 2 



— 1— / 
T l ds 



2s 



E 



> 0. 



(32) 



where Z is the geodesic deviation field (Jacobi) orthogonal to the geodesic 
one. When the point-universe bounces against the potential walls, it is 
reflected from a geodesic to another one thus making each of them unstable. 
By itself, the positive Lyapunov number ( |32"| ) is not enough to ensure the 
system chaoticity, since its derivation remains valid for any Bianchi type 
model; the crucial point is that for the Mixmaster (type VIII and IX) the 
potential walls reduce the configuration space to a compact region (Th), 
ensuring that the positivity of A„ implies a real chaotic behaviour, provided 
the factorized coordinate transformation in the configuration space 



a 



a (9,0 



P+ = e/M& + (0,0 , /?_ = e /M&_ (6,0 , (33) 



where /, a, b± denote generic functional forms of the variables r, 8, ^. 



4 Statistical Mechanics Approach 

In order to reformulate the description of the Mixmaster stochasticity in 
terms of the Statistical Mechanics, we adopt in (|l4] ) the restricted time 
gauge t' = 1, leading to the variational principle 

sj (Pt§+Pe§- u adm) df = 0. (34) 
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In spite of this restriction, for any assigned time variable r (i.e. n) there ex- 
ists a corresponding function / (t) (i.e. a set of Misner-Chitre-like variables 
able to provide the scheme presented in Section ^) defined by the (invertible) 
relation 

# = ZW 7V( ) 2/ ( } 

dr 12D v ; v ' 

As a consequence of the variational principle ([$3) we have again the expres- 
sion ©. 

In agreement with this scheme, in the region Tjj where the potential 
vanishes, we have by (|35|) cITIadm / df = 0, i.e. e = E = const, (by other 
words the ADM Hamiltonian approaches an integral of motion) . 

Hence the analysis to derive the invariant measure for the system follows 
the same lines presented in || ||. 

Indeed we got again a representation of the Mixmaster dynamics in terms 
of a two-dimensional point-universe moving within closed potential walls and 
over a negative curved surface (the Lobachevsky plane ||]), described by the 
line element (j3ll). Due to the bounces against the potential walls and to 
the instability of the geodesic flow on such a plane, the dynamics acquires 
a stochastic feature. This system, admitting an "energy-like" constant of 
motion (e = E), is well-described by a microcanonical ensemble, whose Li- 
ouville invariant measure reads 

dg = A5 (E — e) d^dOdp^dpQ , A = const. (36) 

where 5 (x) denotes the Dirac function. After the natural positions 



= -y|f=y cos </> , p e = e\Ji 2 - lsin^, (37) 

being < < 2 7r, and the integration over all possible values of e 
(depending on the initial conditions, they do not contain any information 
about the system chaoticity), we arrive to the uniform invariant measure 

du = d£d9d(f>-^ . (38) 

The validity of our potential approximation is legitimated by imple- 
menting the Landau-Raichoudhury theorem near the initial singularity 
(placed by convention in T = 0, where T denotes the synchronous time, 

5 Such a theorem, within the mathematical assumptions founding Einsteinian dynam- 
ics, states that in a synchronous reference it always exists a given instant of time in 
correspondence to which the metric determinant vanishes monotonically. 



s 



i.e. dT = —N (r) dr), we easily get that 23 vanishes monotonically (i.e. for 
T — ► we have dlnD/dT > 0). In terms of the adopted time variable r 
(23 — > => /(r) — > oo), we have 



din 23 dlnDdT 



dT 



dT dT 



din 23 
dT 



iV(r) 



(39) 



and therefore 23 vanishes monotonically for increasing r as soon as, by (|3J 
dr/dr > 0. 

Now the key point of our analysis is that any stationary solution of the 



Liouville theorem, like (15), remains valid for any choice of the time variable 
r; indeed in || the construction of the Liouville theorem with respect to the 
variables (£, 8, <p) shows the existence of such properties even for the invariant 
measure (f3S|). 

We conclude by remarking how, when approaching the singularity / —* oo 
(i.e. TCadm E), the time gauge relation (pq) simplifies to 



df 

dr 



Ee -2f+^eS 

12 



N (t) e 



-2/. 



(40) 



in agreement with the analysis presented in ||, during a free geodesic mo- 
tion the asymptotic functions £ (/) , 6 (/) , <p (/), are provided by the simple 
system 



d£ 
df 



£ 2 — 1 COS i 



dO 
df 



sm< 



Once getting £ (/) as the parametric solution 
P 



df 



£ sin< 

7F^ 



(41) 



sin ( 



pcos </> arctanh ' 1 



sin </> 2 +a 2 (l+cos 2 ) 
2 ap cos </> 



p = \l a 2 + sin 2 « 



ap cos i 



a,b = const. £ 3i 



+ 6 



(42) 



it reduces, for a free geodesic motion, equation ( f40| ) to a simple differential 
one for the function / (r). 

However, the global behaviour of £ along the whole geodesic flow, is de- 
scribed by the invariant measure ( |38| ) and therefore relation fl40| ) takes a 
stochastic character. If we assign one of the two functions / (r) or N (/) 
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analytically, the other one acquires a stochastic behaviour. We see how the 
one-to-one correspondence between any lapse function N (77) and the asso- 
ciated set of Misner-Chite-like variables, ensures the covariant nature with 
respect to the time-gauge of the Mixmaster universe stochastic behaviour. 

We are very grateful to Remo Ruffini for his valuable comments on this 
subject. 
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